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ABSTRACT 

This  report  is  a  compendium  of  properties  of  a  mathematical 
function  (the  5  function)  which  arises  in  the  modeling  of  turbulent 
boundary- layer  flows  near  solid  walls.  The  e  function  is  formally 
defined  as  a  triple  integral  which  cannot  be  easily  expressed  in 
terms  of  elementary  or  known  special  functions.  Thus  series  and 
asymptotic  expansions  are  developed  here  as  well  as  a  table  of  inte¬ 
grals  (most  of  which  are  indefinite).  The  5  function  and  its  deriva¬ 
tive  are  tabulated  correct  to  ten  significant  figures.  A  FORTRAN 
function  routine  for  the  calculation  of  e(x)  and  ='(x)  is  also 
given. 
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1.  Introduction 


The  diffusion  equation 


au  _  a2u 

at  ~  ay7  * 


(1.1) 


arises  in  a  wide  variety  of  applications  in  engineering  including 
fluid  mechanics  and  heat  and  mass  transfer.  In  many  situations,  the 
relevant  solution  of  equation  (1.1)  is  a  function  of  the  similarity 
variable. 


=  _X_ 


(1.2) 


When  the  solution  of  equation  (1.1)  is  required  for  y  (or  n)  in 
the  semi-infinite  range  0<y<°°,  a  boundary  condition  is  normally 
known  at  y=0  as  well  as  some  particular  type  of  asymptotic  behavior 
as  y-*».  Specific  types  of  asymptotic  conditions  which  can  occur 
might  include  the  following  situations: 

i)  u  must  approach  zero  or  a  constant  exponentially,  viz. 

u  -  A0  ,  as  y  +  -  ,  (1.3) 

where  A0  is  a  constant  (or  zero); 

i i )  u  must  approach  zero  or  a  constant  algebraically,  viz. 

u  ~  A-|  +  A2f(t)rf01  as  n  00  ,  (1.4) 

where  f(t)  is  a  function  of  t,  A-]  and  are  con_ 
stants  and  a>0;  and 

iii)  u  must  become  large  algebraically  for  large  y,  viz. 

u  -  A3g(t)n6  +  ...  as  n  ,  (1.5) 

where  g(t)  is  a  function  of  t,  A^  is  a  constant  and 
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For  the  types  of  asymptotic  conditions  described  by  equations  (1.3), 
(1.4)  and  (1.5),  the  solution  to  equation  (1.1)  can  often  be  written 
in  terms  of  parabolic  cylinder  functions  or  alternatively  in  terms 
of  repeated  integrals  of  the  error  function  (Abramowitz  and  Stegun, 


i 

i 


i 

t 

i 

( 


1964). 


A  rather  different  type  of  asymptotic  condition  for  large  y 
arises  in  the  study  of  velocity  and  temperature  profiles  in  the  near- 
wall  region  of  turbulent  boundary-layer  flows  (Black,  1968;  Scharnhorst, 
Walker  and  Abbott,  1977;  Scharnhorst,  1978;  Weigand,  1978;  Walker, 
Scharnhorst  and  Weigand,  1986);  here  a  solution  of  the  diffusion  equa¬ 
tion  is  required  whose  time-average  over  a  finite  interval  of  time 
is  logarithmic  for  large  y.  It  is  in  this  context  that  the  ~  func¬ 
tion  arises. 

The  details  of  the  turbulence  model  discussed  by  Walker, 

Scharnhorst  and  Weigand  (1986)  are  complex  and  will  not  be  discussed 
here.  However  an  essential  feature  of  the  analysis  is  that  a  solu¬ 
tion  of  equation  (1.1)  is  required  subject  to  the  asymptotic  boun¬ 
dary  condition. 


u  -  B  In  y  +  C  as  y  -*■  °°  ,  (1.6) 

where  B  and  C  are  known  constants.  The  solution  of  equation  (1.1) 
satisfying, 

u(0,t)  =  A,  (1.7) 


where  A  is  constant,  is  given  by. 


U  =  [C  +  2)  -  V  +  £rt(.t)}]  erf(n)  +  Aerfc(n) 

+  —  H  (n)  .  (1.8) 

Here  y  is  Euler's  constant  and  the  5  function  is  formally  defined 

as  the  triple  integral, 

n  5  5 

*  2  2  f  .  2 

n(n)  =  e”^  e?  e_t  dtdcd?  .  (1.9) 

J0  0  jo 

The  function  defined  by  equation  (1.9)  cannot  be  easily  expressed 
in  tern  of  elementary  or  known  special  functions.  Consequently 
because  the  function  appears  in  velocity  and  temperature  profile  rep¬ 
resentations  of  the  flow  in  the  near  wall  region  of  turbulent  boundary 
layers  (Walker,  Scharnhorst  and  Weigand,  1986),  it  is  worthwhile  to 
develop  a  compendium  of  useful  properties  of  this  function  and  this 
is  the  purpose  of  this  report. 

The  plan  of  the  report  is  as  follows.  In  section  2,  a  number  of 
general  properties  of  the  5  function  are  given  and  h(x)  and  its  deriva¬ 
tive  s'(x)  are  plotted  in  figures  1  and  2  respectively.  The  notation 
is  generally  consistent  with  Abramowitz  and  Stegun  (1964)  and  is 
defined  in  a  separate  section  on  page  iv;  generally  x  is  used  to 
denote  a  real  variable  and  a  denotes  a  positive  real  constant.  In 
section  3,  a  number  of  indefinite  and  definite  integrals  are  listed; 
for  the  indefinite  integrals,  the  constants  of  integration  have  been 
omitted.  In  three  of  the  integrals  given  in  section  3,  a  new  function 
arises  which  has  been  denoted  by  erfei(x);  some  properties  of  this 


function  are  given  in  Appendix  A.  Finally  a  FORTRAN  program  for 
the  evaluation  of  e  and  e'  is  given  in  Appendix  B  as  well  as  tabu¬ 
lated  values  of  these  functions  correct  to  ten  significant  figures 
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2.  GENERAL  PROPERTIES 


2.1  Definitions 


fx  ,2  f5  2  rz  _t2 
S (x)  =  I  e  ^  e  dtdcd?, 

J0  Jo  Jo 

/—  (*x  r2  2 

- (x)  =  -y  e"^  e^(c)dcd£, 


2  fx  2  rc  t2 
E'(x)  =  e  x  e  1  dtdc, 

J0  0 


2  x  2 

’  (x)  =  5  e"x  ec  eA^(c)dc. 
£  n 


2.2  Differential  Equation 


E"  +2x5'  =  $ 


j-eAi(x),  5(0)  =  E '  (0)  =  0. 


2.3  Integral  Representation 


sw  =  --ij?p  /4- 


*5“*  Un(p)  +  £]dp 


2.4  Symmetry  Relations 


iWS) 


e"x2  "  2Jdl_  X2J+1 


2.5  Series  Expansions 
e"x2 

E(X)  =  4  (2 j  +  1 ) ! i 

Here  d .  =  y  +  iHj  +  1 )  =  f  k"1 . 

J  k=l 


-x2  » 

_  e  r  2W 


H'(X)  =  4  j(2j2-  1")!!  x2j 

co  (-1)J+1c.  2i+] 

5(X)  =  2  (2j  +  l)j!  x 


H*(x)=l  l  H^Ux* 


j  =  l 


J!  J 


X I  <  *. 


X,  <  °°1 


X|  <  °°, 


xj  < 


In 


v  1  J 

(2.11)  and  (2.12),  c.  =  in  (2)  +  \  |  <Mj  +  1/2)  =  \  (2i-l) 

J  2  2  .=1 


(2.9) 


(2.10) 

(2.11) 

(2.12) 


-1 


2.6  Asymptotic  expansions 

H(x)  ^  X  f  +  ?  "  \  l 

L  j=l 

(2j  -  1 ) !  !"| 

X  -*■  °D 

(2.13) 

- 1  /  \  „  »/tF  v  (2j  -  1 ) ! !  ^ 

i(x)^ 4  jlo  ■ 

x  ■+  co. 

(2.14) 

2.7  Expression  as  a  series  of  inconplete  gamma  functions 


* 

! 


-  (x)  = 


V  TT 


8  j=i  jr(j  +  1/2) 


r\  Y  ( J  +  1/2  . x"). 


(2.15) 


¥ 


for  negative  argument  the  double  factorial  is  defined  to  be  identically 
equal  to  one. 


2.8  Relation  to  Dawson's  Integral 


- '  (x)  =  y-  eA|$(x)F(x) 


2.9  Maximum  of  E 1 (x) 


='(1.379005630....)  =  0.3048918974 


5  (1.379005630....)  =  0.2257293172 


yji*. 


3.  INTEGRALS 


3.1  Combination  of  H  and  E'  with  Powers 


1.  E(ax)dx  =  xE(ax)  +  ^  E'(ax)  -  'j-  xe^(ax)  -  ^  e 


.2  2 

-a  x 


» 

!.  xE 


E(ax)dx  = 


<2a2*V  1}  5(ax)  +  K-H'(ax) 
4aa  ,a 


2  2 
x  -a  x 

16a  e 


(2a2x2-l  )eA{(ax) 


3.  xnE(ax)dx  =  p,(a,x)S(ax)  +  p2(a,x)E 1 (ax)  +  p3(a,x)e4|$(ax) 


p4(a,x)e 


J  2 
-a  x 


In  3^\_  3  for  n  a  positive  even  integer,  viz.  n  =  2k  with  k  =  0,1,2,. 


n  _  x _ 

P1  "  2k+l 


L 


»  Pp  2k+l  I  -i  I  (ax) 

2  2 (2k  +  l)a2k  j=0  J- 


✓?  k!  f  a  J 
4(2k  +  l)a2k  j=0  (2<j  *  ^ 


a2j~  2j+l 

7 TTTT  a 


p4  -  - 


- k!  ,kt1  I  i  (ax)2j(c.+.  c . ) 

4(2k  +  l)a2k  1  j=0  J'  kl  J 


where  c.  is  defined  in  (2.12)  for  j  >  1  and  cn  =  0.  In  3.1  3  for  n  a  positive  S 

j  0  v 

odd  integer,  viz.  n  =  2k  +  1  with  k  =  0,1 ,2,3, . . . :  S 


Y3W* 


3.2  Combinations  of  5  and  with  Exponentials  and  Powers 


x,  <  oo. 


Here 


.  (-nJt1  f  mc 

j  j!  • 


where  is  defined  in  (2.12). 


2. 


roo  2 

e~x  H‘(ax)dx  =  i  [ -J —  £n(a2  +  1) 
0  8  7  a  +1 


2  _  x 

3.  xe”x  5 1  (x)dx  “  -  |  H * (x)  +  j|-  vif  (x). 


f  xVx  S'(x)dx  -  -  |  e'x  ='(x)  -  ^|e‘x  c*J(x)  +  e*<(/EO 


e"x  5 '  (x)dx. 


x3e”x  S'(x)dx  =  -  ^  (2x2  +  l)e"x  E'(x)  +  ca,i'?(x) 


/if  -x2  ,/  »  1  -2x2 

-  yg-  xe  e'lfHx)  -  e 


32 


6. 


xne"x  5 '  (x  )dx  =  -  jj-  e"x  5'(x)  +  y-  xn‘1e~x  e4)((x)dx 


+  x"-2e-x  5.(x)dx< 


For  the  evaluation  of  the  second  term  on  the  right  side  of  equation  3.2 
6  the  following  reduction  formulae  ore  useful: 


-13- 


Iv‘v\V 


xn~lg_x  ^(xjdx  ~  -  j  e  x  eAjJ(x)  +  JL  xn”2e~2x  dx 


+  lniil 
2 


xn-3e"X  e/L^(x}dx, 


•  x"-V2*2dx  •  -  f'e-2*2  +  x"-3e-2x2dx> 


^  -t2_ 


e"c  H(t)dt  =  l  b.x2j+2, 

j=l  J 


where 


4TJ+T7T  1i, 


and  c.  is  defined  in  (2.12). 


I.  |Vt2S(t)dt  =  |-  3j£w(2), 


where  G  is  Catalan's  constant. 


9. 


-x2  e~x  1  -x2 

xe  x  5(x)dx  ~  j  H(x)  +  j  e  x  S' 


(x)dx. 


10. 


n  -x2  xn_1  -x2  1  n-1  -x2 

x  e  x  5(x)dx  -  -  j  e  X  H(x)  +  j  x  e  x  5‘(x)dx 


+  indl 
2 


xn"2e"x  5(x)dx. 


For  the  evaluation  of  the  second  term  on  the  right  side  of  \2_  10  see 
3.2  6. 


-14- 


where 


for  n  > 


where  b 


rawrow  r 


.  is  defined  in  3.2  1. 


L  2 

12.  j  e  t  ='  (at)  y-  =  yg  -  an.c.ot{a.)\j\  -  aAcot( a)],  -»  <  a  < 


13.  (V^H'tt)  %  *  - 
'0  t 


i  ='<*>- 4 


0  8  * 


♦5 


rx 


.  ^2(t) 

J0  t: 


14. 


e  1  :'(t)  y  s  y  ln{  1  +  /?)  -  £u ( 2 ) , 

0  t  4 


15. 


r°°  _t2  j  •. 
e  r  H 1  (t)  —  = 


-x 


r°°  4.2 


(n-l)x 


CT5'(x)  -  7STTT 


e-^  =< [+.\ 

x  -  U)  t"’2 


^  r.-t2 


dt 


+  2lCTFj  * 

X  t 


16. 


x  4.2 

e 

0 


t2=(t)  dt  =  y  iiitll  b  x2j+l 
’  t  2j+l  Djx 


Xi  < 


.  is  defined  in  3.2  7. 


17. 


0<x>  2 

e_t  H(t)  ^  =  0.047128  94918. 


for  n  >  2. 


3.3  Combinations  of  5  and  5'  with  the  Error  function  and  Powers 

2 


1. 


(x)E 1  (x)dx  =  CA|$(x)E(x)  - 


e"x  S(x)dx, 


2. 


» 

2 

-X 

eA.j{(x)E(x)dx  = 

XCA^(x)  +  ~ 

H(x)  +  J  c^(x)E'  (x) 


2_ 

Sit 


j  e_t  5' (t)dt 

^  [W(x)  +  ^  e"x  eA6(x)  e^(/2x)J 
3.  xeA)J(x)E '  (x)dx  =  -  j  eAj$(x)E'(x)  + 


_t2 

e  L  E ' (t)dt 


+  xeA^2(x)  +  i  e’x  cai$(x)  -  '-j-  e^(/2x) 


4. 


xe^(x)~(x)dx  = 


+  ^xe"x 


H(x) 


+  ^  XCA($(x)  +  j- 


-X 


H’(x)  -  y|  XW(X) 


1  -x2  ,  v  1  -2x2 

a  xe  eAj((x)  -  —7=  e 


5.  j  x"^(x)='  (x)dx  =  xneAj$(x)H(x)  -  JL 


x  e"x  E(x)dx 


-  n 


xn"^eA($(x)5(x)dx. 


For  evaluation  of  the  second  term  on  the  right  side  of  3.3  5  see  3.2  10 


3.4  Combinations  of  E  and  E1  with  Logarithms  and  Powers 


rx 


1. 


Cn(t)E'(t)dt  =  £n(x)E{x)  -  f  ^4^-  dt, 
0  1 


rx 


£n(t)E(t)dt  = 


2(x)  +  i='(x)  -  £xe/i<(x>  -  !■'*  +  } 


-X 


[cn(x)  -  1 J  +  j-  xe4j$(x)  +  1  je‘x  -  lj  -  1  E^x2) 


-  4  ln(x)  ~  q  ~  2 


rx 


'(^T- 


3. 


t£n(t)E ' (t)dt  =  - 


^  E'(x)  -  —  xe*$(x)  -  \  e"x  +| 


ln() 


+  h 


H'(t)  f-  ^xexj(x)  -i(e'x  -  1) 


+  £  [E1 (x2)  +  2 in(x)  +  y], 
rx  P  2 

Un(t)S(t)dt  =  S(x)  +  x  E 1  (x)  -  5  (2x2  -  l)eA|fl 


4. 


3*5  Miscellaneous  Integrals 


1. 


2. 


V^eA^atJdt  =  e3^ 1  (ax) , 

0  ^ 


x  2  2 

e*  eAj$2(t)dt  =  ^  ex  eA^(x)='(x)  -  £  E(x), 


fX  t2  ,3 


2  _x  __  ,2, 


Qe  ^  (t^dt  /jf  e  ^6  (x)='(x)  -  -  e/Lrf(x)S(x) 

+ ^  re't2=(‘)dt. 


4. 


f*  t2  ,  2 


e1"  2A)$c‘’(t)eAj5(t)dt  =  JL  E '  (x)  +  ^  ca^c  (x)E(x) 


16 

7r/rr 


x  *2 

e_t  E(t)dt, 


fx  t2  o  2 

e  eA($c(/2t)eA(${t)dt  =  -j=  zsifc  (/2x)ex  E'(x) 


V? 


x  _t2 

e  1  E ' (t)dt. 
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APPENDIX  A:  THE  FUNCTION  ERFEI(x) 

In  integrals  3.1  10,  3.1  ll  and  4,  the  function  defined  by 
equation  (A.l)  arises.  In  this  Appendix,  some  useful  properties  of  this 
function  are  given;  the  function  is  plotted  in  figure  A.l  and  tabulated 
in  table  A.l 

A.l  Definition 


e**ec(x)  =  fX  dt. 

n  1 


A. 2  Integral  Representation 


1  ~2.2 


O  1 

zA^eJ-(x)  =  -  —  x  e  £n(t)dt, 

/7r  Jn 


X  I  <  ». 


A. 3  Symmetry  Relation 


eA&<U(-x)  =  -  eA^(x) 


A. 4  Series  Expansion 


e*6ex.(x)  =  y:  l  — — 2 —  x* 
^  j=0  (2j  +  1)  j! 


X  <  oo. 


A. 5  Asymptotic  Expansions 


<lk^qx{x)  'V  ln{ 2x)  +  £  +  — 


2  23  22  1689 

x2  4?  j?  16x8 


4881  .  1 

^  -J’ 


t 


K*. 


H 


A. 6  Indefinite  Integrals 


1. 


eA^ex(x) dx  =  xeA^(x)  -  xeA^(x) 


2. 


xneA)$ex.(x)dx  =  xn+1eA^eX(x) 


1 


(n+1  )/n 


Y(1  + 


TABLE  A. 1 


X 

eAfiziix) 

X 

eA^zi{x) 

.00 

0.00000  00000 

2.00 

1.67535  04327 

.10 

0.11271  27665 

2.10 

1.72395  62837 

.20 

0.22468  00025 

2.20 

1.77036  53206 

.30 

0.33518  26238 

2.30 

1.81475  13483 

.40 

0.44355  25900 

2.40 

1.85727  26708 

.50 

0.54919  39999 

2.50 

1.89807  27340 

.60 

0.65159  93946 

2.60 

1.93728  10996 

.70 

0.75036  04607 

2.70 

1.97501  45953 

.80 

0.84517  28573 

2.80 

2.01137  85246 

.90 

0.93583  54106 

2.90 

2.04646  78598 

.00 

1.02224  43601 

3.00 

2.08036  83701 

.10 

1.10438  36667 

3.10 

2.11315  76561 

.20 

1.18231  25702 

3.20 

2.14490  60827 

.30 

1.25615  16257 

3.30 

2.17567  76063 

.40 

1.32606  83591 

3.40 

2.20553  05035 

.50 

1.39226  35013 

3.50 

2.23451  80088 

.60 

1.45495  85267 

3.60 

2.26268  88709 

.70 

1.51438  49639 

3.70 

2.29008  78382 

.80 

1.57077  57076 

3.80 

2.31675  60822 

.90 

1.62435  83548 

3.90 

2.34273  15672 

4.00 

2.36804  93746 

Values  of  vifa t(x)  for  x  >  4  (correct  to  ten  significant  figures)  may 
be  obtained  by  using  the  first  two  terms  of  A. 5. 


APPENDIX  B:  FORTRAN  FUNCTION  ROUTINE  AND  TABULATED  VALUES  OF  5(x)  AND  H'(x) 

In  Table  B.l,  the  listing  of  a  function  subprogram  is  given  for  the 
numerical  evaluation  of  S(x)  and  5'(x).  The  functions  5(x)  and  5'(x)  are 
accessed  from  the  calling  program  by, 

Y  =  XI(X,EPSI ,IERR)  ,  Y  =  XIP(X,EPSI  ,IERR) , 
respectively.  In  these  calling  statements 

X  =  the  argument  of  either  H(x)  or  E'(x), 

EPSI  =  the  assigned  tolerance  which  is  dependent  of  the  number 
of  significant  figures  desired  for  E(x)  or  E'(x), 


IERR  =  an  integer  variable  returned  by  the  function  subprogram. 

(IERR  =  0  if  the  tolerance  EPSI  is  met;  IERR  =  1  other¬ 
wise.) 

The  subprogram  accepts  only  zero  or  positive  values  of  x  and  returns 
a  zero  value  if  x  <  0.  If  the  evaluation  of  5(x)  or  H ' (x )  is  required  for 
x  <  0,  the  calling  statement  should  reflect  the  appropriate  symmetry  relation 
(either  (2.7)  or  (2.8)). 

For  x  <  5.0  the  routine  evaluates  E(x)  and  E'(x)  using  the  power  series 
expansions  (2.9)  and  (2.11)  in  the  DO  loops  commencing  with  cards  20  and  29, 
respectively.  The  power  series  in  (2.9)  and  (2.10)  are  summed  until  they 
have  converged  to  the  number  of  significant  figures  dictated  by  the  assigned 
value  of  the  tolerance  EPSI;  once  convergence  occurs  control  is  transferred 
out  of  the  DO  loops  to  statements  3  or  5,  respectively.  The  DO  loops  commencing 


with  cards  20  and  29  are  initially  set  (by  card  11)  so  that  at  most  100  terms 
in  either  power  series  will  be  evaluated.  If  convergence  occurs  before  the 
DO  loops  are  satisfied,  a  reliable  result  has  been  obtained  and  I  ERR  =  0  is 
returned  to  the  calling  program.  If  the  DO  loops  are  satisfied,  the  returned 
value  of  - (x)  or  5'(x)  may  not  be  reliable  to  the  specified  tolerance  (EPSI) 
and  IERR  =  1  is  returned  to  the  calling  program. 

For  x  >_  5.0,  the  asymptotic  expansions  (2.13)  and  (2.14)  are  used  for 
the  evaluation  of  5  and  The  sums  in  (2.13)  and  (2.14)  are  evaluated  in 
the  DO  loops  commencing  with  cards  38  and  46,  respectivley.  Control  is  trans¬ 
ferred  out  of  these  DO  loops  when  the  last  term  added  to  these  sums  is  less 
than  the  desired  exit  tolerance  EPSI;  again  a  returned  value  of  IERR  =  1 
indicates  this  condition  has  not  been  met  after  100  terms  in  the  series  have 
been  summed. 

Computed  values  of  5  and  are  given  in  Table  B.2;  these  values  were 
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evaluated  with  EPSI  =  10  and  are  correct  to  ten  significant  figures  (the 
last  digit  is  rounded  off).  The  CPU  time  required  to  produce  the  entire 
table  was  0.842  seconds  on  a  CDC  6500.  The  storage  required  for  XI  is  171 
words.  The  program  was  written  for  a  CDC  machine;  for  computers  which 
operate  with  less  presision,  XI  may  easily  be  modified  for  double  precision 
if  a  high  degree  of  accuracy  is  needed. 


■ 

1 

FUNCTION  XI<X,EPSI«IERR) 

1 

.V 

OAt*  SRPI.GAMO/l, 772453850905516, 0.57721566490 153286/ 

2 

XI=OslF=l 

3 

l>. 

IF ( X , LE • 0  .  )  RETURN 

4 

.■4’ 

60  TO  1 

5 

•X 

entry  xip 

6 

J 

XI=0.STF=2 

7 

IF(X.LF.O.)  return 

8 

■<; 

X2=X*X 

9 

FAC=2.*X2 

10 

m=ioo 

11 

ICRRsO 

12 

SUM=Q, 

13 

>2 

SUMT=0. 

14 

I 

TERMsi. 

15 

>*, 

IF(X.ge.5.)  60  TO  110 

16 

IF ( IF .EG, 2 )  60  TO  100 

17 

TERHsX 

18 

alpha=i. 

19 

./> 

jt 

00  2  1=1, M 

20 

term=term*fac/float<2*i*i) 

21 

K 

SUMsSUK+TERM*  alpha 

22 

IF ( ABS ( ( SUM-SUMT ) /SUP )  ,LT ,EPSI }  60  TO 

3 

23 

ALPHAsALPHA+1. /FLOAT ( 1*1 ) 

24 

”  »  ! 

SUMTsSUM 

25 

IERR=1 

26 

■ . 

XI=0.25*EXP<-X2I*SUM 

27 

• 

return 

28 

•  mJ 

00  4  1=1, M 

29 

TERM=TERM*FaC/FL0AT<2*I-11 

30 

$0 

SUM=SUM«-TERM/FL0aT  < I ) 

31 

Ss 

IF <  ABS ( ( SUM- SUM T ) /SUN ) ,LT • EPS  I )  60  TO 

5 

32 

SUMTsSUM 

33 

lERRcl 

34 

XI=0,25*EXP(-X2) *SUM 

35 

Return 

36 

i 

IF<IF.EQ.2>  60  TO  120 

37 

DO  6  1  =  1, M 

38 

TERM=TERM*FL0ATI2*I-1)/FAC 

39 

TER«A=TERM/FLOAT(I) 

40 

/,■] 
»  • 

IF ( TtRMA.LT , EPSI )  GO  TO  7 

41 

w  m  \ 

.  «* 
0  a.  1 

sum=sum*terma 

42 

K 

IERR=1 

43 

• 

XIsSRPI»(  ALOGIX2)*GAMO-SUM>/6. 

44 

RETURN 

45 

00  8  1=1, M 

46 

:V«i 

TERMsTEKM*FLOAT( 2*1-1) /FAC 

47 

■*.> 

IFITER4.LT. EPSI )  GO  TO  9 

4  6 

o 

SUM= SUM* TERM 

49 

IERR=1 

50 

M 

-  ■  * 

XI*SRPI*U.*SUM)/I4.*X) 

51 

return 

52 

ENO 

53 

•'A 

TABLE  B.2 


X 

(x) 

—  1 

(x) 

X 

- 

(x) 

0.00 

0.00000 

00000 

0.00000 

00000 

5.00 

0.83647 

55948 

0.10 

[-3)0.16600 

18215 

(-2)0.49667 

94096 

5.10 

0.84543 

37558 

0.20 

[-2)0.13122 

31678 

(-1)0.19474 

75580 

5.20 

0.85421 

05888 

0.30 

[-2)0.43419 

16955 

(-1)0.42390 

90215 

5.30 

0.86281 

34482 

0.40 

[-1)0.10012 

91790 

(-1)0.71967 

91837 

5.40 

0.87124 

92399 

0.50 

[-1)0.18885 

34957 

0.10603 

37513 

5.50 

0.87952 

44570 

0.60 

[-1)0.31290 

22308 

0.14221 

74031 

5.60 

0.88764 

52166 

0.70 

,-1)0.47321 

19062 

0.17818 

00344 

5.70 

0.89561 

72816 

0.80 

[-1)0.66848 

42292 

0.21182 

20008 

5.80 

0.90344 

60950 

0.90 

[-1)0.89550 

98227 

0.24144 

25483 

5.90 

0.91113 

68005 

1.00 

0.11496 

21681 

0.26583 

80743 

6.00 

0.91869 

42648 

1.10 

0.14252 

14171 

0.28433 

51804 

6.10 

0.92612 

30976 

1.20 

0.17162 

64540 

0.29676 

42616 

6.20 

0.93342 

76696 

1.30 

0.20168 

03441 

0.30338 

65749 

6.30 

0.94061 

21293 

1.40 

0.23212 

96272 

0.30479 

17529 

6.40 

0.94768 

04176 

1.50 

0.26249 

14660 

0.30178 

34527 

6.50 

0.95463 

62820 

1.60 

0.29236 

94322 

0.29526 

87699 

6.60 

0.96148 

32891 

1.70 

0.32145 

89291 

0.28616 

31134 

6.70 

0.96822 

48360 

1.80 

0.34954 

41840 

0.27531 

74919 

6.80 

0.97486 

41616 

1.90 

0.37648 

92054 

0.26347 

05669 

6.90 

0.98140 

43556 

2.00 

0.40222 

51381 

0.25122 

40720 

7.00 

0.98784 

83683 

2.10 

0.42673 

61718 

0.23903 

76697 

7.10 

0.99419 

90188 

2.20 

0.45004 

56933 

0.22723 

80552 

7.20 

1.0004 

59002 

2.30 

0.47220 

38365 

0.21603 

69560 

7.30 

1.0066 

30898 

2.40 

0.49327 

70769 

0.20555 

32997 

7.40 

1.012Z 

17176 

2.50 

0.51334 

00923 

0.19583 

58849 

7.50 

1.0187 

20202 

2.60 

0.53246 

98095 

0.18688 

40873 

7.60 

1.0246 

42244 

2.70 

0.55074 

13614 

0.17866 

52362 

7.70 

1.0304 

85478 

2.80 

0.56822 

55924 

0.17112 

81812 

7.80 

1.0362 

51993 

2.90 

0.58498 

77346 

0.16421 

31835 

7.90 

1.0419 

43795 

3.00 

0.60108 

69077 

0.15785 

86265 

8.00 

1.0475 

62811 

3.10 

0.61657 

61587 

0.15200 

51934 

8.10 

1.0531 

10895 

3.20 

0.63150 

28204 

0.14659 

81709 

8.20 

1.0585 

89830 

3.30 

0.64590 

90301 

0.14158 

84599 

8.30 

1.0640 

01331 

3.40 

0.65983 

23067 

0.13693 

27628 

8.40 

1.0693 

47051 

3.50 

0.67330 

61180 

0.13259 

32920 

8.50 

1.0746 

28582 

3.60 

0.68636 

04075 

0.12853 

72415 

8.60 

1.0798 

47458 

3.70 

0.69902 

20615 

0.12473 

61723 

8.70 

1.0850 

05160 

3.80 

0.71131 

53151 

0.12116 

54032 

8.80 

1.0901 

03113 

3.90 

0.72326 

20998 

0.11780 

34521 

8.90 

1.0951 

42696 

4.00 

0.73488 

23379 

0.11463 

15474 

9.00 

1.1001 

25240 

4.10 

0.74619 

41927 

0.11163 

32117 

9.10 

1.1050 

52029 

4.20 

0.75721 

42797 

0.10879 

39140 

9.20 

1.1099 

24305 

4.30 

0.76795 

78461 

0.10610 

07801 

9.30 

1.1147 

43268 

4.40 

0.77843 

89246 

0.10354 

23j34 

9.40 

1.1195 

10081 

4.50 

0.78867 

04646 

0.10110 

83963 

9.50 

1.1242 

25865 

4.60 

0.79866 

44457  1 

[-1)  0.98789 

72500 

9.60 

1.1288 

91708 

4.70 

0.80843 

19760  1 

[-1)  0.96578 

07181 

9.70 

1.1335 

08663 

4.80 

0.81798 

33776  1 

[-1)  0.94465 

96875 

9.80 

1.1380 

77749 

4.90 

0.82732 

82618  I 

[-1)  0.92446 

64968 

9.90 

1.1425 

99954 

- 1  (x) 

(-1)0.90513  96674 
(-1  )0.88662  31905 
(-1)0.86086  59127 
(-1)0.85182  10068 
(-1)0.83544  55111 
(-1)0.81969  99290 
(-1)0.80454  78777 
(-1)0.78995  57812 
(-1)0.77589  25981 
(-1)0.76232  95830 

(-1)0.74924  00732 
(-1)0.73659  93004 
(-1)0.72438  42224 
(-1)0.71257  33729 
(-1)0.70114  67264 
(-1)0.69008  55781 
(-1)0.67937  24349 
(-1)0.66899  09175 
(-1)0.65892  56716 
(-1)0.64916  22887 

(-1)0.63968  72325 
(-1)0.63048  77738 
(-1)0.62155  19304 
(-1)0.61286  84121 
(-1)0.60442  65712 
(-1  )  0.59621  63571 
(-1)0.58822  82742 
(-1)0.58045  33441 
(-1)0.57288  30702 
(-1)0.56550  94057 

(-1)0.55832  47240 
(-1)0.55132  17913 
(-1)0.54449  37413 
(-1) 0.53783  40521 
(-1) 0.53133  65247 
(-1) 0.52499  52632 
(-1) 0.51830  46560 
(-1  )  0.51275  93591 
(-1  )  0.50685  42801 
(-1  )  0.50108  45635 
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